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Abstract. The pulsing inertial oscillation (PIO) model is a nonlinear, time-dependent, 
translating vortex solution of the inviscid, compressible fluid dynamic equations in the middle 
troposphere. The translation of this vortex during a pulse is strikingly similar to that of a 
supercell storm — a rotating thunderstorm that can generate tornadoes and hail. Two studies 
were performed to test the hypothesis that some supercell storms are manifestations of a PIO 
pulse. The first study applied the model to an intense interior draft whose buoyancy was 
bounded by a temperature excess of ± 12 K . The peak updraft speed achieved was 41.5 m s' 1 
and the peak Rossby number was 92.9. The study also pointed to an advanced concept for 
attaining higher values. The second study applied the PIO model to a supercell storm as a whole 
and succeeded in replicating its bulk properties, such as mesocyclonic circulation, net mass and 
moisture influxes, and time track. This study also identified a critical feature of the PIO model 
that could be tested against storm data: The average vertical draft is downward before the turn in 
the storm track and upward afterwards. In the conventional theory, the average vertical draft is 
upward from storm inception until dissipation. These differing draft predictions were compared 
with the best available data, which are surface mesonetwork data. These data were found to 
support the PIO model. However, surface data alone are not conclusive, and further 
measurements are warranted. 
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1. Introduction 

A supercell storm is a rotating thunderstorm that can generate tornadoes and hail. A salient 
feature of supercell storms is that, upon spin-up, the track of a cyclonic supercell turns to the 
right and that of an anticyclonic supercell turns to the left (Fujita and Grandoso [1]). We present 
a novel solution for a tropospheric vortex that has the same feature. The model is nonlinear and 
is called the pulsing inertial oscillation (PIO). Two studies are made to test the hypothesis that 
some supercell storms are manifestations of a PIO pulse. 

These studies revive Ferrel's [2] concept that the rotation of a supercell storm is caused 
by the Coriolis force. This concept has remained dormant for many years, partly because of 
scaling studies by Wallace and Hobbs [3], Morton [4], Holton [5, pp. 1-19], Davies-Jones [6], 
and others that concluded that the Coriolis force is too weak to spin up a supercell storm in the 
observed time. The PIO model is a counterexample that negates this conclusion. 

The source of rotation in the conventional theory, as reviewed by Bengtsson and Lighthill 
[7], Kessler [8], Ray [9], Klemp [10], Church et al. [11], and others, is horizontal vorticity 
associated with vertical wind shear. An initial buoyant updraft tilts and stretches the horizontal 
vortex lines and forms two counter-rotating storms that split apart, as shown numerically by 
Schlesinger [12]. Although such storm splitting has been observed on radar by Fujita and 
Grandoso [1], Brown [13], and others, it occurs infrequently compared with cyclonic supercells, 
according to Davies-Jones [6]. This apparent difficulty with the conventional model was 
resolved as follows: Maddox [14] showed that the vertical windshear vector in the lower 
troposphere usually turns clockwise with increasing height in the vicinity of a supercell storm. 
Using numerical simulations, Klemp and Wilhelmson [15] then showed that this turning vector 
strengthens the cyclonic storm and weakens the anticyclonic storm. If the weakened storm 
should produce little precipitation, it would be invisible to radar. 

Although the conventional model is highly developed and successful in certain respects, 
some fundamental questions remain unanswered. For example, as mentioned by Lilly [16] and 
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Klemp [10], the model does not explain the rightward translation of a cyclonic supercell. Also, 
the conventional model misplaces the initial development of a tornadic vortex - a vortex of 
Rossby number about 1000 that precedes a tornado by about 20 min. According to the 
conventional model, a tornadic vortex originates at the surface, while Doppler radar indicates 
that it originates at midheight in a supercell [6,11,17]. 

The PIO model consists of a partly analytic, partly numerical solution of the inviscid, 
compressible fluid dynamic equations for a translating vortex in the middle troposphere. 
Compatible solutions are assumed for the boundary and upper layers. The core radius of the PIO 
model is arbitrary below about 60 km, so the model can be applied to a supercell storm as a 
whole or to an intense internal draft. The latter is the subject of our first study, which is aimed at 
determining the upper limits to updraft speed and vorticity relative to the production of hail and a 
tornadic vortex. 

The second study treats the storm as a whole. In its present rudimentary form, the PIO 
model does not resolve coexisting updrafts and downdrafts, but instead treats the time-dependent 
spatial average vertical draft. This average vertical draft is found to be negative during spin-up - 
- in contrast to the conventional model, where it is positive during spin-up. We compare these 
differing model predictions with the best available data, which are found to be surface 
mesonetwork data. 

Costen and Stock [18] applied the PIO model to a supercell hailstorm that occurred 
during the National Hail Research Experiment (NHRE) in northeast Colorado, USA, on 22 July 
1972. (We will refer to this supercell, which is documented by Foote and Fankhauser [19], as 
the Colorado hailstorm.) The mesocyclone was already spun up before this storm entered the 
NHRE mesonetwork, so no analysis of surface data was done. 

In the second study, we apply the PIO model to a supercell hailstorm that originated 
inside the mesonetwork of the Cooperative Convective Precipitation Experiment (CCOPE) in 
southeast Montana, USA, at 0030 UTC on 12 July 1981. (We will refer to this storm, which is 
documented by Fankhauser et al. [20], as the Montana hailstorm.) The CCOPE mesonetwork is 
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described by Knight [21] and the data quality by Wade [22]. After originating within the 
mesonetwork, this storm developed into a cyclonic supercell that produced baseball- to 
grapefruit- size hail, and its track turned to the right while still on the mesonetwork. Therefore, 
this data set is appropriate for testing the PIO and conventional models. 

2. PIO model 

2.1. T ANGENTIAL COORDINATE FRAME 

The coordinate frame is a local tangential Cartesian frame where the origin is fixed at mean sea 
level (MSL) and the x—, y—, and z - axes point eastward, northward, and upward, respectively. 
(Cylindrical coordinates r.0. and z are also used, where x = rcos6 and y = r sin0.) 
Although the x - and y - axes do not curve with the surface, this frame is adequate for 
describing tropospheric flows within a horizontal radius of about 60 km [23, pp. 164-166]. This 
radius is sufficient to contain the convective core of a supercell, but not all the outer flow. 

In the Cartesian frame, the fluid velocity is denoted by v = (u,v,w) (in m s'l), the 
vorticityby co = curl v = (^, 77 , £) (in s' 1), the divergence by D = div v 
= (du/d x + dv/dy + dw/dz ) (in s'l), and the angular velocity of the Earth by 
O = (0,£2_ v ,£2,) (in rad s' 1). The components O v and Q , although slowly varying functions 

of y , are treated as constants. 

2.2. S TORM IDEALIZATIONS AND ZEROTH-ORDER APPROXIMATION 

Miller et al. [24] present Doppler radar measurements of a supercell storm that show a boundary 
layer near the surface, a highly sheared layer near the top, an updraft that has its maximum at 
midheight, two downdrafts, and a mesocyclone. As mentioned, the cylinder shown in Figure 1 is 
intended to represent, in highly simplified form, either the entire convective core of a supercell 
or an intense internal draft. Our analysis is confined to the midtropospheric layer (b < z < h), 

which excludes the boundary and upper layers wherein compatible flows are assumed. Lateral 
entrainment, friction, heat conduction, and radiation are neglected. 

As shown in Figure 1, the PIO model is a rotating and dilatating vertical draft of radius 
a(t ) . The centerline is located at (x r ( t ), y c (t)) . The mesocyclone is represented by a Rankine 
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vortex of vorticity £(t) , which is uniform inside the core (r < a(t)) and zero in the outer 
region ( r>a(t )). To emphasize compressibility, the density p (in kg nr 3 ) is taken to 
decrease exponentially with height z , which also implies that the temperature lapse is zero 
between b and h . For our purposes, this isothermal approximation is preferable to the 
constant density approximation used in the supercell studies by Rotunno and Klemp [25] and 
Lilly [26]. Rising parcels of air expand with divergence D(t ) , which is also uniform inside the 

core and zero in the outer region. The core fluid is uniformly buoyant with a normalized density 
deficit B(t ) , which is called simply the buoyancy. The fields co(t), D(t ) , and B(t) are zeroth- 

order approximations to the instantaneous spatial average vorticity, divergence, and buoyancy in 
the core. This property will be used when the model is compared with storm data. 

The density of the outer fluid at radius a is given by 

p°(a,z) = p i V £Z (la) 

where the fields in the outer region are distinguished by the superscript 0 and Z = z~b . The 
coefficient p b ° and the inverse density scale height £ (in nr ' ) arc constants. Inside the core, 

the density is given by 


p(z,t) = p h °[l-B(t)]e- eZ (lb) 

The buoyancy B(t) « 1 , because it corresponds to a temperature excess or deficit of at most 
about 12 K [24] and terms of 0(B 2 ) will always be neglected. 

The continuity equation is given by 


— + \Vp + pD =0 

c )t 


( 2 ) 
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Substituting the inner density (lb) in (2) and retaining the effects of B(t) on the buoyant force 
but neglecting its effects on D(t ) , in accordance with one of the Boussinesq approximations 
[27], we have 

D(t) = e w(t) (3a) 


It follows that the vertical component of velocity w(t) is also spatially uniform inside the core 
and that D(t) is purely horizontal divergence given by 


, du dv 
D(t) = — + — 
dx dy 


(3b) 


Equation (3a) is a key approximation of this study and w(t) and D(t)/e will be used 
interchangeably. Also included in the core is a horizontal flow with vertical windshear that is 
represented by the constants U b and V h and the horizontal vorticity components %(t) and 
ri(t). 


3. Inner reduction 

3.1. Velocity field 

The complete velocity field inside the core is given by 

u(x, y,z,t) = U„ +rj(t)Z + ^\D(t)X(x,t)-C(t)Y(y,t)] (4a) 


v(x, y, z, t) = V b -Z(t)Z + 1 [. D(t)Y(y , t) + £(t)X(x, t)] 


(4b) 


w(x,y,z,t) = w{t ) = 


Pit) 

£ 


(4c) 


where X(x,t) = x-x c (t) and Y(y,t) = y - y c (t) . The constants U b and V h represent the 
stationary part of the horizontal flow inside the core, and i](t )Z and -^(t )Z represent the 
inner vertical windshear. The terms - C ( t)Y / 2 and £(t)X / 2 represent the inner flow of the 
Rankine vortex. The angular velocity of the core fluid 0 (in rad s" ' ) is given by 
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9 = S- 


where an overhead dot denotes an ordinary time derivative. The terms D(t)X / 2 and D(t)Y / 2 
represent the inner radial horizontal flow that results from the divergence. 

3.2. Momentum equation 

The inviscid momentum equation is given by 


dv „ v 2 V p „ 

— + mxv + V — + -^- + V<T + 2nxv = 0 (6) 

dt 2 p 

where v 2 is v ■ v , p is the pressure (Pa), is the geopotential, gz (m^ s"2), g is the 

gravitational acceleration (9.81 m s'2). 

Solving (6) for Vp and substituting the inner velocity field (4) we find 


|( = - £(G,-G,v + Du) 

(7a) 

f =- P (<;,+Gj, + Dv) 
dy 2 

(7b) 

dp ( D 

3z e J 

(7c) 

where 


G,=DX -CY -Dx +C y c +2riZ + 2 — (ri + 2Q. v ) 

£ 

(7d) 

G,=Dr + gX-(i'-Dy,- 2i;Z-2-t; 

£ 

(7e) 

G 3 

(7f) 


3.3. Second-order partial derivatives and nonlinear harmonic equations 

Pressure jumps are known to occur at gust fronts in the surface layer (z < b) , as described by 

Haltiner and Martin [23, pp. 305-307]. However, in the midtropo spheric layer, we assume that 
p and its first- and second-order partial derivatives are continuous functions of v, y, and z 
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(except at radius a, where p must be continuous but its derivatives can be discontinuous). It 
follows that for (r < a,b < z ^ h ) 

d 2 p _ d 2 p d 2 p _ d 2 p 

ckdy dydx ’ d.xdz dzd.x 


d 2 i 


d 2 ; 


dydz dzdy 


( 8 ) 


Substituting the cross derivatives of (7a), (7b), and (7c) in (8), dividing by p , and setting the 


coefficients of X , F, and Z 

individually to zero in each equation, 

we obtain the following set 

of coupled ODE's: 

| = i|n(f+4a,)-z>S] 

(9a) 


+ 4£1,)+Di;] 

(9b) 


i=-D(g + 2aJ 

(9c) 


D = i[c(f+4Qj-C 3 ] 

(9d) 


- ^ ( y, - v b )+ ^ fc - u J- n - - a , - + ia z v„ = o (9e) 

Ui c -U,)+ky c -V„)+l;--n,^-2n : U b =0 (9f) 

2 2 £ £ 


Equations (9e) and (9f) can be solved for the translational velocity of the centerline 
(in m s' 1 ) 


= u„ + - a„ + — 2 - [£(-£ w + 2 a u „)+ D(ri w - 2 ay b )] 


y rl 


£ 2 +D 2 


(10a) 
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y c ~Vb + ^2 +d 2 [c( _7 7 w + 20 _V b )+£>(- £ w + 20 z U b )] (10b) 

At the equator where O v is maximal, the term 20. v / £ in (10a) contributes about 1 m s'l to 

the translation. 

3.4. Pressure field 

To obtain the inner pressure field p , we first substitute (9) and (lb) in the pressure gradient (7), 
which gives 


^ = -p,V' z (l-fl)£L- (lla) 

OX £ 

^=p„V' z (l-B)£2,T (lib) 

oy £ 


dp 

dz 


-pb ) e~' z {l~B)\g + w-20 s 


U b +r,Z + ^(DX-CY) 


(He) 


These partial derivatives can now be integrated to obtain p as 


p(x,y,z,t) = p b e (l B) 8 \ 1 + 

^ + U b +riZ + UDX-^Y)f 

«l « 8 

|_£ 2 JJ 


( 12 ) 


This result shows that the PIO model has no pressure drop at its center, where X = Y = 0 . A 
2-mb pressure drop was measured in the Montana hailstorm. According to Davies- Jones [6], the 
pressure deficit in mesocyclones is typically about 5 mb. Such a pressure deficit indicates a 
deviation from the PIO model. Reasons for such deviations are discussed in Section 5.4. 
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3.5. Balance of forces 

Equations (9e) and (9f) can be multiplied by pA , where A = k a 2 , and written in vector form to 
give a balance of forces per unit length (in N m' 1 ) 

z x (r. -Y b )+^AD( r c -V b )~ pAcoxz w-AV H p-2pA£2xzw-2pAQ z zxV b =0 

(13a) 

where carats denote unit vectors, r c = x x c + y y c , \ h = xU h + yV h , and, by (11a), (lib) and 
(lb) 


V H p = —xClp — + yClp — = zxClp — — zx(zxCl)p — (13b) 

£ £ £ £ 

The first term is analogous to the Magnus force or Kutta-Joukowski force on a cylinder 
with a bound vortex of circulation AC, ; differences are a factor of 1/2 and the definition of V,, 

as the stationary part of the fluid velocity inside the core. The second term is analogous to the 
force on a source tube of strength AD . The third term is the force exerted by the inner 
horizontal vorticity on the vertical draft. The fourth term is the force exerted on the core by the 
horizontal pressure gradient. The fifth term is the Coriolis force on the vertical draft. The last 
term is the Coriolis force on the stationary part of the flow inside the core. 

Thus far, we have reduced the inner solution to a time-dependent set of ODE's, except 
that we have no equation for B(t) or a(t ) . These quantities are determined by the jump 

conditions at the cylindrical interface between the inner and outer solutions, as shown in Figure 
1. The cross section is circular only if the model does not translate. For simplicity, we will 
determine Bit) and a(l) for a nontranslating model and assume that these results apply 

approximately to the translating model. This assumption will fail if the inner flow does not 
block the relative environmental wind (Brandes [28]), which will then erode the core. 
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4. Nontranslating PIO model 

4.1. Inner solution 

For the PIO model to remain centered on the origin (/.<?., for x c = y c = 0), we must take 
U h =V h = q = 1 ] = Q . , = 0 . The velocity (4) and pressure (12) then become in cylindrical 

coordinates 

v r (r,t) = ^p-r (14a) 


v g (r,t) 



(14b) 


v z (t) = w(t) = 


m 

£ 


(14c) 


p(z,t) = p b °e~ eZ -[l-5(0] 
£ 



w(t) 

8 


(14d) 


4.2. Outer solution 

With the Rankine vortex as a model, we take the external vorticity ox' and divergence D° to 
be zero. Hence, the outer velocity field is given by (for r > a) 


0/ ^ D(t)a\t) 
(r,t) = 


(15a) 


V e (. r,t ) 


gt)a 2 (t ) 

2 r 


(15b) 


v : ° = w° =0 


(15c) 
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The field v r is continuous at the interface r = a(l ) , as required by the jump conditions for the 
continuity and momentum equations. The field v e is also continuous, but w, p , D, and C, all 

have finite discontinuities. The requirement that the interface move with the fluid gives 



(16) 


The core buoyancy Bit) can be determined from the jump condition that the pressure be 
continuous at r = a. To determine the outer pressure field, we first solve the momentum 
equation (6) for the pressure gradient in cylindrical coordinates and then substitute the outer 
velocity field from (15) and the tendencies from (9c), (9d), and (16) to obtain 

dp° _ § F P {) 

dr £ r 


z, \ 

1 a 

2 

r 


(17a) 


dp 0 

"do 


= 0 


(17b) 


dp 0 
dz 



(17c) 


where F , the Froude number, is given by 


(17d) 


As with the inner solution, we require 


2 0 

d n 


-\2 0 

a p 


drdz dzdr 


(r > a,b < z ^ h ) 


(18) 
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Consequently, the outer density p° must satisfy 


5p! 

dr 



£ r y 



(19) 


The solution of this equation that remains bounded at infinity and satisfies the boundary 
condition (la) at r = a is given by 


P°(r,z,t) = p b 


M 

F 

^ F 

r 


exp 

— £ Z H 

2 

‘“I 


( 20 ) 


Substituting this result in (17) and integrating for the outer pressure we have the hydrostatic 
result 


p°(r,z,t) = —p ° (r, z, t) = p b ° — 


f a ) 

F 

^ F 

r 


exp 

— £ Z H 

2 

M 


( 21 ) 


Because of the smallness of F (Table 1), p° and p° are essentially independent of r within 
the 60 km radius of validity. As r->°o,p° and p° slowly approach zero. 

4.3. Buoyancy and inner pressure 

Equating the outer pressure (21) and the inner pressure (14d) at r = a we obtain 


5(0 = ^ 

8 

This expression is similar to that given by Darkow [29] for the vertical acceleration of a 
nonentraining buoyant parcel. The inner pressure is now independent of t and given by 


P(z) = p h 


0 —e~ eZ 

£ 


(23) 


By (9d), the buoyancy (22) can also be written 
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B = ^-[^+AQ. z )-D 2 } (24) 

2 eg 

4.4. Supplementary formulas 

Although the model is now complete, formulas for certain additional quantities will be useful for 
comparison with meteorological data. The mesocyclonic circulation T is given by (in m2 s' 1) 

T = n a 2 £ (25) 

The net mass influx M at cloud base b is given by (in kg s' 1) 

M =7ta 2 wp h °{l-B) (26) 

As mentioned earlier, the outer temperature T° in the model is a constant; its value is 
given by 

t ° = JL ( 27 ) 

eR 


where R, the gas constant, equals 287 J kg 1 K'l. The inner temperature T is given by 

T(t) = T°[l + B(t)] (28) 

Thus, the model core is also isothermal, but its temperature varies with the buoyancy. 

The energy equation, as given by Holton [5, p. 51], can be written 

fa 

q = pc {j, 

where q (in W m~3) is the diabatic heating rate, which in the PIO model is due only to latent 
energy release, and c p equals 1004 J kg' 1 K 1 . Substituting the inner pressure (23) and 

temperature (28) and integrating over the core, we obtain the total heating rate Q (in W) 
required to support the PIO 


+ v • V T — P v ■ V p 

dt 
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Q = P, 


o 71 Cl' 


c p T B + w 


dT 

c ^ +s 


[\-e £(hh) ] 


(30a) 


Although the model troposphere is isothermal between b and h , we have retained dT / dz to 

be more accurate when applying this formula to an actual supercell. The corresponding net 
water vapor influx M v (in g s'l) is given by 


M. 


Q 

4 


(30b) 


where L c , the latent heat of condensation, equals 2.5 xlO 3 J g 1. 

5. Nondimensionalization and analytic solution 

5.1. Nondimensionalization 

We define dimensionless (primed) variables by 


t'= 20. j 


(31a) 




(31b) 


(where we will use C,' interchangeably with the Rossby number Ro ) 


B'- 


egB _£g{T-T°) 


20 


2 0 2 T° 


a = 


(31c) 
(3 Id) 


where subscript 0 denotes an initial value at t'= 0. 

The model Equations (9a)-(9d), (24), (16), and (5) can then be written in dimensionless 


form, as follows 
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2^=n'(C'+2 )-£>'§■ 

dt 

(32a) 

2 C ^ 7 = -£'(C+2)-D‘ri' 
dt' 

(32b) 

^ = — D'(C'+ 1) 

dt' 

(32c) 

2^ = C'(C'+2)-D' 2 
dt' 

(32d) 

B'=£'{£'+2)-D' 2 

(32e) 

da' D’a' 
dt' ~ 2 

(32f) 

dd _ C 
dt' ~ 2 

(32g) 


5.2. Analytic solution 

Equations (32c) and (32d) uncouple from the other equations. Although this pair is nonlinear, 
the equation for dD' / d C is linear in D' 2 and leads to an analytic solution. Subject to the 
initial conditions D'( 0) = 0 and C'(0) = Co' > the analytic solution is given by 


C '«') 


cost'- 


iol 

Co’+2 


1 + 


fo'(fo'+2) 


■cos t 


(33a) 


D\n 


Co '(Co '+2) Sint' 

I Co '(Co ’+2) sin t'l 


(Co- CO 


C r . A 

C'+-^- 

Co'+l 


V 


(33b) 
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Equations (32c), (32d), and (32f) yield, with a'(0) = 1 


a\f) 


ICq'I 

7c +D a 


(33c) 


from which it follows that the Froude Number F , as given by (17d), is a constant of the motion. 
By (32a), (32b), and (32f), we have ^' 2 +ri' 2 = -^yj^ 0 ' 2 +r] 0 ' 2 , which shows that the magnitude 

of the inner horizontal vorticity vector is inversely proportional to the core radius. Integration of 
(32g) gives, with 0(0) = 0 


0(0 


Co ' +1 

Co'+l 


- arctan 

Co '+1 tan 

|7Y 

- 

_ 

l 2 )\ 


2 


(33d) 


Plots of Z)' versus C and D\ B\ a', and 0 versus t' are shown in Figure 2 for 
C 10 ' = -0.814 (Montana hailstorm simulation) and C 20 ' = -1.186. In the (C',D') phase plane, 
there are equilibrium points at (0,0) and (-2,0) . The point (0,0) is the center for an oscillation 
that gives cyclonic pulses and (-2,0) is the center for anticyc Ionic pulses. The line £'= -1 

separates these two oscillations and represents infinite contraction and spin-up for both. Since 
initial values are chosen along the line D'= 0 , two values of £ 0 ' give the same solution, 

except that the pulses are displaced by n. (Had we chosen £ 10 '= 4.39 and £ 20 '= - 6.39 , the 
same pulses in Figure 2 would have occurred at t'= 0, ± 2k,... ) Cyclonic and anticyc Ionic spin- 
up both occur in a contracting downdraft. The downdraft reverses in the middle of a pulse and 
spin-down occurs in an expanding updraft. This draft reversal results from a pulse of positive 
buoyancy. 

5.3. Extrema formulas , etc. 

Values of the fields at t'=0andt'=7T are related as follows: 
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C'(7T) 

Co ' = 


_ -Co' 
Co'+l 

-CW 
C'OO + i 


(34a) 

(34b) 


D 0 '=D\n) = 0 (34c) 

Z? 0 '=Co’(Co'+2) (34d) 


BW=z ISr (34e) 

fl’(^) = |Co'+l| (340 


The fields D', B' , and 0 cyclomc have extrema that occur at intermediate times: 


D' 


max 


c;(c+ 2) 

2(f„'+l) 


(34g) 


D' 


B' 


-D' 


8(fo'+l) 2 


(34h) 

(34i) 


e „ = 


arctan^o'+l - arctan — = 


V c o'+i 


(£ 0 '> -l) cyclonic pulse 


0 min = - 0 max (C o ' > - 1 ) cyclonic pulse 


(34j) 

(34k) 


As the intensity of a cyclonic pulse increases, 0 max — > 7T / 2 ; hence the total angular excursion of 
the core fluid (0 max - 9 mm ) -> n . 

The dimensionless time interval to spin up a cyclonic peak from £' = 0 is given by 


A t's P in- up = arccos 


f r y \ 

bo 

Co’+2 


(osf;) 


(341) 
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As the pulse amplitude increases, the spin-up time decreases. This feature is typical of nonlinear 
oscillations [30]. 

Since the limiting factor on the amplitude of a pulse is the core buoyancy, it is useful to 
rewrite some of the preceding relations in terms of the initial buoyancy B 0 ' , where B 0 '> — 1, as 

follows: 


1 + tJi + B 0 ' (cyclonic pulse) 

(35a) 

l-y]l + B 0 ' (anticyclonic pulse) 

(35b) 

B\k)= ~ B ° 

1 + B 0 ' 

(35c) 

\ b A 

D' - L^J 

“ 2^1 + *„' 

(35d) 

_-(2 +B lt y 

min 8(1 + B 0 ') 

(35e) 


Equation (35e) shows that as the pulse amplitude gets large, the positive peaks of B'(t') become 
8-times greater than its negative peaks. 

5.4. Maximum vorticity and updraft speed 

In this section, we apply the PIO model to an intense vertical draft in the interior of a supercell 
storm. The objective is to see if the PIO can reproduce a hail-producing updraft of w ~ 50 m s ' 1 
or a tornadic vortex of Rossby number Ro ~ 1000 . As mentioned, we limit the temperature 
excess in the core {t -T 0 ) to the range (-12K<r-r° < 12 k). The parameters for the 
Montana hailstorm are listed in Table 1: 0_=5.2623x 10" ^ s' 1 , £= 1.19x 10"4 m' 1 , and T°= 287 
K. Taking T -T° = 12 K in (31c), we obtain 5' max = 5 0 '= 8.81xl0 3 . Equations (35a) and 
(35d) then give Co' = C'max = R° m ax = 92.9 and D' max = 46.9 , from which we obtain 
Cmax = 9.78 x 10” 3 s' 1 and w max = 41.5 ms' 1 . 
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This value for w max is 17 percent short of the goal. The value for £ max , although triple 
the maximum vorticity of 3 x 1 0 3 s' 1 measured in the interior of the Montana hailstorm (Section 
7.1), is an order of magnitude less than the goal of a tornadic vortex. However, the 8 to 1 
asymmetry between positive and negative peaks of buoyancy in result (35e) has indicated 
something new to try: In the preceding calculation, we took the minimum dimensionless 
buoyancy B' min = -1. lx 10 3 , for which T -T° = -1 .5 K . If, instead, we take 
B ' m jn = -8.81x 10 3 , which corresponds to the negative limit T — T° = -12 K , then the positive 
extremum 5' max will be far off the allowable scale. Since the PIO goes through negative values 
of buoyancy first, it could proceed until it encounters the positive limit (T - T" = 12 K) . Then 
the PIO would stall, leaving the internal draft in a state of prolonged downflow, contraction and 
spin-up - a scenario that should produce a higher Rossby number. A different stall scenario 
could increase the updraft speed and also explain the 5 mb central pressure deficit typically 
observed for mesocyclones. After any stall, however, the problem requires a different 
formulation that is beyond the present scope. 

6. Physical nature of PIO 
6.1. Inertial flow 

We present a simple physical picture for the two oscillations shown in Figure 2. The model is 
nontranslating. The period n / Q , (in s), which is half the period of a Foucault pendulum, 

suggests that the oscillations are inertial. As described by Holton [5, pp. 64-65] and Haltiner and 

Martin [23, pp. 178-181], inertial flow is frictionless horizontal flow in which a fluid parcel 
moves with constant speed v, along an arc such that the centrifugal force balances the Coriolis 

force. The horizontal pressure gradient V H p is zero - which agrees with the absence of a 
central pressure deficit in the PIO model. The radius r t of the arc is given by (in m) 
r t = v, I2Q and the motion is anticyclonic. Since we take Q constant, the parcel traverses a 
circle, as shown in Figure 3(a). Introducing the dimensionless coordinates r '= r t / a 0 and 
v ; '=v,. /2Q z a 0 we have r,'= v, ' . 
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6.2. Core contraction and spin-up 

We now consider the fields in Figure 2 at t'= 0 . Figure (3b) is an aerial view of the 

dimensionless core at this time. The velocity of parcels on the periphery is given by 

v, ' = — £ 0 7 2 . The inertial radius of these parcels, r f '= — £ 0 7 2 , is also shown. The physical 

explanation for the cyclonic pulses in Figure 2 is as follows: Figure 3(c) shows the inertial arcs 
for 4 To '= — 0.8 14 during the first half period (0 < t'< k) when the core contracts and spins up 

cyclonically. As the core contracts, the parcels are compressed by descending in a negatively 
buoyant downdraft. The negative buoyancy results from evaporative cooling of moisture 
deposited in the middle troposphere by a previous updraft. This downdraft reverses at t'=n 
because of a pulse of positive buoyancy from condensational heating of moist low level air that 
is lifted by the gust front of the downdraft. Figure 3(d) shows the complete cycle including spin- 
down in an expanding updraft. Although we have only followed parcels on the periphery, every 
parcel in the core simultaneously traverses a similar inertial circle. 

The anticyclonic pulses in Figure 2 have a similar explanation. Figure 3(e) shows the 
inertial arcs for £ 0 '= -1.186 , where the core contracts and spins up anticyclonically. The 

vertical speed and buoyancy are the same as for cyclonic spin-up. A singularity occurs as 
rj'—> 0.5 (or -l) that corresponds to infinite contraction and spin-up. In the 

troposphere, this mathematical singularity is avoided because of physical limitations on the 
buoyancy. 

6.3. Blocking meso- anticyclone 

Finally, if r t '= 1 (or £ 0 '= -2), we obtain a steady meso-anticyclone. (The full set of model 
equations shows that x c = 2X2 v / £ , y c =D = w=B = Q = 0, and that a,%,r),U b , and V b are 

arbitrary constants, subject to the constraint that the core not be eroded. The meso-anticyclone 
translates eastward at about one-tenth the speed of weather patterns and, therefore, would have a 
blocking effect.) 
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7. Application to Montana hailstorm 

7.1. Storm description 

Our second study compares the PIO model with data from the Montana hailstorm. This storm 
initially appeared on radar inside the CCOPE mesonetwork at about 0030 UTC on 12 July 1981. 
The midpoint of the right turn in its track occurred at 0200 UTC. Between 0206 and 0209 UTC, 
three Doppler radars measured an interior vorticity maximum of 3x 10' 3 s' 1 . At 021 1 UTC, the 
supercell was producing baseball- to grapefruit- size hail. At about 0250 UTC, it moved off the 
mesonetwork and evolved into a mesoscale convective complex. We hypothesize that the PIO 
started sometime between 0030 and 0100 UTC, proceeded through one pulse, and stopped 
shortly after 0250 UTC. 

7.2. Initial values and parameters 

The initial values and other parameters used in the PIO model to simulate this hailstorm are 
listed in Table 1, where t = 0 coincides with the midpoint of a pulse. Dimensionless plots for 
D' , B' , a', and 9 are already given by the cyclonic pulses in Figure 2 with a trivial shift of 

the t' -origin. Besides these plots from the analytic solution, we numerically integrate (9a), (9b), 
and (10) to obtain and x c (t), y c (t ) , which are plotted in Figure 4. (The numerical 

procedure is a fourth order Runge-Kutta-Fehlberg code. The time step is automatically adjusted 
within the range 10' ^ to 1, subject to an error tolerance of 10"6 per step.) The five solid 
triangles in Figures 2, 3(d), and 4 relate the simulation to universal time during the hailstorm. 
The average vertical draft in the model is also tested qualitatively against surface mesonetwork 
data at these five times. 

The initial values in Table 1 were chosen as follows: As shown in Figures 2(b) and 4(b), 

the midpulse is coincident with the center of the turn and, therefore, occurs at 0200 UTC. The 
minimum core radius a Q was chosen to include the 40 dBZ low level radar echo at 0200 UTC 

so that the core would capture the major convection as well as the mesocyclonic vorticity. The 
initial value £ 0 , which determines the degree of nonlinearity, was chosen so that a(l) 

replicates the growth rate of the radar echoes after 0200 UTC, as given by [20]. The initial 
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values Co > Vo and the constants U b , V h were chosen to simulate the time track of the 
mesocyclone. The PIO time track from Figure 4(b) is shown superimposed on the observed time 
track in Figure 5, where fairly good agreement is evident. 

With these initial values, the model computes the peak mesocyclonic circulation T as 
2.8 x 105 m^ s' * at 0200 UTC. Although T was not reported for the Montana storm, the 
computed value compares with a typical value of 5 x 10^ s"' , as given by Davies-Jones [6]. 

The net mass influx M at cloud base is computed as 1.32 x 10^ kg s' 1 at 0230. This quantity 
was not reported either; however, the computed value compares well with a measured value of 

1.3 x 1 09 kg s'l for M during a similar time in the Colorado hailstorm. The model computes 
the net water vapor influx M v at 0230 as 4.2 x 10^ g s'l. This value compares with a 

measured value of 8.5 x 10^ g s'l — again from the Colorado hailstorm. According to our 
simulation, M v continues to rise after 0230. We speculate that the PIO stops shortly after 0250 

UTC partly because the net water vapor influx required to maintain it becomes too large. 

7.3. Comments on time track, etc. 

The eastward drift of the kidney-shaped track shown in Figure 4(b) results from the term 
2 Q y I £ in (10a). The upper (westward) track segments that occur between the pulses are not 

observed for supercells and are included only to show the complete oscillation. These upper 
segments suggest a second reason for the PIO to stop shortly after a pulse. As the track turns 
toward the north and eventually toward the west, the relative environmental wind increases. The 
inner flow can no longer block the relative environmental wind, which then erodes the core. 
Incidentally, we take h = 6.3 km in Table 1 because a vertical radar section appears to show 
substantial core erosion above this height [20]. 

We seek an explanation for the right turn in the track at 0200 UTC. Figure 4(a) shows 
that the average horizontal vorticity vector rotates clockwise about *4 turn from 0130 to 0230 
UTC with changing from negative to positive at about 0200 UTC. The average vertical draft 

w also changes from negative to positive at 0200, as shown in Figure 2(c). Equation (10b) for 
y c confirms that these sign reversals in w and g cause the right turn in the track. (If the 
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model is run with =rj Q = 0 , no turn at all occurs at 0200 UTC.) Newton and Newton [31], 

Rotunno and Klemp [25], and others have also linked the right turn in the track to horizontal 
vorticity but not to a reversal of the average vertical draft. 

Besides the time track, Figure 5 shows the model core at the five times. Figure 6 shows 
the corresponding low level radar echoes. The core approximately matches the 40 dBZ radar 
echo at 0200 UTC. However, at earlier times, 0130, 0100, and 0030 UTC, the model core is 
substantially larger than the corresponding radar echo. We conclude that much of the average 
downdraft predicted by the PIO model during this earlier period would have been invisible to 
radar. 

(The individual radar scans shown in Figure 6 were begun at 0030:36, 0102:05, 0132:05, 
0202:19, and 0232:14 UTC on 12 July 1981. The scan at about 0030 UTC was taken at 0.5° 
elevation by the NCAR CP-2 radar (10-cm wavelength) located at (0,0) km on the x,y coordinate 
frame shown. (The origin of the coordinate frame is located at 46° 25' 56" N lat and 105° 56' 
10" W long.) The remaining scans were taken at 1° elevation by the Bureau of Reclamation 
Skywater radar (5-cm wavelength) located at (2. 5, -2. 2) km. The reflectivity grayscale shown to 
the right of each figure goes from 10 to 60 dBZ.) 

8. Comparison of the PIO and conventional models with surface mesonetwork data 

8.1. Critical data 

The salient feature of the PIO model is the contracting average downdraft that spins up the 
mesocyclone with help from the Coriolis force. This downdraft occurs before the track turns to 
the right — before 0200 UTC for the Montana hailstorm. After 0200, the average draft becomes 
upward. In the conventional model, the average vertical draft is upward from storm inception 
until dissipation. 

Relevant data before the right turn are not easy to find. If the downdraft should contain 
little condensed water, it would be invisible to the eye and camera as well as to radar. This 
downdraft could be detected by aircraft wind measurements; but, as shown by Foote and 
Fankhauser [19], it takes aircraft about 45 min to circumnavigate a storm. Once a probable 
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supercell is identified, little time is left for aircraft to reach it and make measurements before the 
right turn. At present, the only way of detecting this downdraft appears to be with a surface 
mesonetwork. The average downdraft would cause an average divergent flow near the surface 
that could be detected by the wind sensors. Conversely, an average updraft would cause an 
average convergent flow. 

Surface data alone are not conclusive. Multilevel wind data between the surface and 
cloud base (1.5 km above the surface for the Montana storm) are needed to measure the net mass 
influx at cloudbase. Such multilevel wind data could be obtained by a mesonetwork composed 
of wind profilers operating in conjunction with Doppler radars. The wind profilers would 
function in the clear air and the Doppler radars in the precipitation regions. Until such data 
become available, however, surface data appear to be the only means for testing the average 
vertical draft sequences in the PIO and conventional models. 

8.2. Unique data set 

After forming inside the CCOPE mesonetwork, the Montana hailstorm traversed the southern 
part of the network where the mesh size was about 20 km and the wind sensors were mounted 6 
m above ground level (AGL). Hence, the term "surface" as we will use it actually denotes 6 m 
AGL. Surface winds were recorded every 5 min at each operating station. The recorded values 
were time averages of data taken during the previous 5 min. Sample vector plots of these 
recorded surface winds were made every half hour. These plots were obtained on 35-mm 
microfilm, entitled "CCOPE PROBE Data — A First Look," from the National Center for 
Atmospheric Research in Boulder, Colorado, USA. 

8.3. Surface data analysis method 

We analyzed these plotted surface wind data at 0030, 0100, 0130, 0200, and 0230 UTC on 12 
July 1981. The analysis was done graphically on 25 by 25-cm prints. A polygon was drawn on 
network points that surround the center of the mesocyclone at each of these five times, as shown 
in Figure 6. The vertices were taken 20 to 40 km from the mesocyclone center. This size was 
needed to accommodate the 20-km mesh size and to capture the model core shown in Figure 5. 
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Some network stations were inoperative or only partially operative when this storm passed 
through; the resulting absence of data caused irregularities in some of the polygons. 

8.4. Average HORIZONTAL divergence versus t 

The average horizontal divergence (Z>) on the surface polygons is plotted versus t in 
Figure 7. Positive (Z>) indicates an average downdraft just above the surface. Similarly, 
negative (D) indicates an average updraft just above the surface. According to Wade [22], 

the error in the surface wind measurements was 0.5 ms'l The error bars in Figure 7 were 
determined by using this value in an error formula given by Davies-Jones [32]. 

Figure 7 supports the PIO hypothesis. The transition from average updraft at 0030 UTC 
to average downdraft at 0100 UTC suggests that the PIO started during this interval. The average 
downdraft continued until about 0200 UTC, when it changed to an average updraft, in qualitative 
agreement with the PIO model 

8.5. Draft structure and time dependence 

Further details of the measured vertical draft were obtained by using interior grid points to divide 
each polygon into triangles. When the horizontal divergence was computed on each triangle, 
coexistent updrafts and downdrafts were resolved. In particular, at 0200 UTC, we found a 
forward- and a rear-flank downdraft. These dual downdrafts are typical of well-developed 
supercell storms, as shown by Miller et al. [24] and Davies-Jones [6]. 

If the measured horizontal divergence on each triangle is multiplied by its area, we obtain 
the horizontal efflux E sfc (in s'l). These effluxes are additive and can be plotted separately 

for the downdraft, updraft, and net draft as done in Figure 8. This figure highlights the 
coexistence of updrafts and downdrafts. The updrafts, subjected to vertical windshear, could 
continuously supply moisture to the middle troposphere for evaporatively cooling the 
downdrafts. Figure 8 shows that the reversal in net draft at about 0200 UTC was caused by a 
rapid growth of the updraft. 
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9. Conclusions 

The pulsing inertial oscillation (PIO) model is a nonlinear, time-dependent, translating vortex 
solution of the inviscid, compressible fluid dynamic equations in the middle troposphere. A 
property of this model is that vertical windshear causes the track of a cyclonic vortex to turn to 
the right during a pulse and that of an anticyc ionic vortex to turn to the left. Such translation is 
typical of supercell storms — rotating thunderstorms that can generate tornadoes and hail. 

Two studies were made to test the hypothesis that some supercell storms are 
manifestations of a PIO pulse. The first study applied the PIO model to an intense internal draft. 
The objective was to simulate a hail producing updraft ~ 50 m s" 1 and a tornadic vortex — a 
vortex aloft of Rossby number ~ 1000 that precedes development of a tornado at the surface. 
The draft buoyancy was limited by practical considerations to a temperature excess of ± 12 K . 
Subject to this limitation, the PIO produced a maximum updraft speed of 41.5 m s' 1 and a 
Rossby number of 92.9. However, the high degree of asymmetry between the positive and 
negative values of buoyancy required to support the PIO suggest that the PIO could easily get 
stalled in a state of prolonged downdraft, contraction, and spin-up. Such a stall could increase 
the Rossby number. A different stall scenario could increase the updraft speed and also explain 
the 5 mb central pressure deficit typically observed for mesocyclones. To explore these 
promising options, however, requires an advanced model that would apply after a stall has 
occurred. 

The second study applied the PIO model to a supercell storm as a whole. The model 
successfully simulated bulk properties such as the mesocyclonic circulation, the net mass and 
moisture influxes, and the time track. Moreover, we found a critical feature of the PIO model to 
test against storm data. The average vertical draft is downward before the turn in the track and 
upward afterwards. This feature of the PIO model contrasts with the conventional model, in 
which the average vertical draft is upward from storm inception until dissipation. 

Surface mesonetwork data were found to be the best data available for testing the vertical 
draft sequences in the PIO and conventional models. Data from a hailstorm that occurred during 
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the 1981 Cooperative Convective Precipitation Experiment (CCOPE) in southeastern Montana, 
USA, support the PIO model. However, surface data alone are not conclusive, and further 
measurements are needed. A mesonetwork consisting of wind profilers operating in conjunction 
with Doppler radars would yield multilevel data that are conclusive. 
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Table 1. Initial values and parameters used in PIO model for 
Montana hailstorm. 
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a Fankhauser el al. [20]. b Above mean sea level (MSL). c Sounding 
from Baker, Montana, at 0110 UTC 12 July 1981 (see also Fig. 5). 
d From (23) at height z = b 
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Figure 1. Pulsing inertial oscillation (PIO) model for convective core of a supercell storm or an 
intense internal draft. 
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(a) Oscillations on the £ - D' phase plane. 

UTC 



(b) Vertical vorticity £ versus t'. 

Figure 2. PIO plots of solutions (33) and (32e) for = ~ 0-814 (Montana hailstorm 
simulation) and ^ 20 = “1-186 with Dq' = 0. Primed quantities represent 
dimensionless core averages; solid triangles relate the simulation to universal time 
during the Montana hailstorm of 12 July 1981. 
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(c) Horizontal divergence D' versus t'. 
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(d) Buoyancy B' versus t'. 


Figure 2. Continued. 
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(e) Core radius a ' versus f. 



(f) Angular displacement of core fluid 0 versus t\ 


Figure 2. Concluded. 


(a) Inertial oscillation of individual air parcel. 



(b) Dimensionless model core at f = 0 and inertial velocity vj' 
and inertial radius of parcels on the periphery. 

Figure 3. Aerial views showing physical nature of PIO oscillations 
in Figure 2. 
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(c) Contraction and cyclonic spin-up for = - 0.814. 



(d) Complete cycle for = ~ 0.814. Solid triangles 
are the same as in Figure 2. 


Figure 3. Continued. 
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(e) Contraction and anticyclonic spin-up for ^q' = —1 .186. 


Figure 3. Concluded. 
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(a) Phase diagram of inner horizontal vorticity components and r\ . 



(b) Time track of model centerline on horizontal plane. 

Figure 4. Numerical solutions of (9a), (9b), and (10) for the Montana 
hailstorm simulation. Initial values and parameters are listed in Table 1 . 
Solid triangles are the same as in Figure 2. 



x, km 

Figure 5. Theoretical and observed time tracks for Montana hailstorm on CCOPE 
mesonetwork. Data are from Fankhauser et al. [20]. Theoretical time track is 
from Figure 4(b), and theoretical core size is from Figure 2(e) and Table 1. 
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(a) 0030 (0030:36) UTC. 



(b) 0100 (0102:05) UTC. 

Figure 6. Low-level radar echoes from Montana hailstorm (12 July 1981) and corresponding 
surface analysis polygons superimposed on CCOPE mesonetwork. The reflectivity grayscale 
goes from 10 to 60 dBZ. Starting times of the radar scans are shown in parentheses. 
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(c) 0130 (0132:05) UTC. 
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(d) 0200 (0202:19) UTC. 


Figure 6. Continued. 
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(e) 0230 (0232:14) UTC. 


Figure 6. Concluded. 
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Figure 7. Average horizontal divergence (D) s j- C of winds measured near 
surface beneath Montana hailstorm versus t. 
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Figure 8. Updraft, downdraft, and net draft components of horizontal 
efflux E s f c computed from winds measured near surface beneath 
Montana hailstorm versus t. 


